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Abstract. - We track the motion of a horizontally vibrated amorphous assembly of bidisperse 
hard disks, for densities ranging across the jamming transition. We derive on very general grounds 
a bound on the dynamical susceptibility in terms of the response of the dynamics to a change in 
density. This generalizes a similar bound recently derived for equilibrium liquids. We find that 
in our experimental system the bound is tight and reproduces the non-monotonic behavior of the 
dynamical susceptibility both in time and density across the jamming transition. The underly- 
ing scaling behavior reveals an intimate connection between anomalous diffusion and dynamical 
heterogeneity. 
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Introduction. — The sudden change from a homoge- 
nous liquid to an amorphous rigid state is a phenomenon 
common to many different systems, among which super- 
cooled liquids, colloids and granular media. This transi- 
tion, called in these different contexts either a glass or a 
jamming transition, has very peculiar properties [1-3]. In 
particular, the associated dramatic increase of relaxation 
times in a very narrow window of temperatures or densities 
is apparently not related to the growth of any static long 
J> range order: simple static correlation functions remain 
quite featureless. On the other hand, growing long range 
I'jjj order has been discovered in dynamical correlations [4,5]. 
A complete characterization of these so-called "dynamical 
heterogenities" has become a central task in the context 
of glass and jamming transitions. Indeed, many theories 
ascribe the slowing down of the dynamics to some kind of 
growing dynamical correlations or cooperativity, starting 
from the seminal work of Adam and Gibbs [6]. Precise 
characterization of the amplitude of these heterogeneities 
has been obtained using the recently introduc;ed dynam- 
ical susceptibility Xi [7,8], mostly in numerical work [4] 
but only in very few experimental situations [9, 10]. Since 
evaluating this susceptibility requires extensive statistics 
on time-resolved density correlators, it turns out to be 
out of reach in most molecular glassy systems. Instead, it 



has been argued that the response of the dynamics to a 
change in a relevant control parameter (density, tempera- 
ture,...) provides the dynamical susceptibility with a lower 
bound [11-13] which potential interest relies on the fact 
that it only involves readily measurable two-point func- 
tions. However, the quahty of this lower bound has never 
been experimentally investigated and remains to be estab- 
lished. Only numerical results on model glass formers are 
available [13]. 

In this letter, we focus on this issue in the general con- 
text of driven athermal systems displaying a jamming 
transition. First, we show that the lower bound to the 
dynamical susceptibility obtained for equilibrium systems 
in [?] can be extended to driven athermal situations un- 
der a very general hypothesis. Second, we provide the first 
experimental test for the quality of this bound in a glassy 
granular monolayer undergoing a jamming transition. We 
directly measure the four-point dynamical susceptibility 
and compare it to its lower bound. We find that this 
bound is tight and reproduces the non-monotonic behav- 
ior of X4 both in time and density across the transition. 

Theoretical context. — Consider a large bulk region 

of size V of a still much larger stationary system. The sys- 
tem can be either in thermal equilibrium or a driven ather- 



p-1 



F. Lechenault et al. 



mal system such as the vibrated grain assembly studied 
below. Wc; assume that the density p = N /V is the only 
conserved quantity, where N is the number of grains in 
the volume V . This is indeed the case in the experimental 
situation we focus on. Because of the local conservation 
law, p evolves on slow time scales, through boundary ef- 
fects. On very long time scales (or on different regions of 
size V), however, the density fluctuates around its mean 
value (p), with fluctuations of order 1/^/V. Take now any 
one-time intensive observable 0{t) and consider the joint 
distribution of 0{t) and p{f), P{0,p) = P{0\p) x P{p). 
In absence of long-range correlations, one expects on gen- 
eral grounds that for large V (i) the probability laws re- 
duce to large deviation functional and (ii) the probability 
law for O in a bulk region far from the boundaries de- 
pend on an external control parameter a only through 
the density, [a denotes e.g. the position of a piston, that 
allows one to change the average density of the system.] 
This can be indeed verified for equilibrium systems and 
for some simple model of athermal systems (see [14] for 
a detailed discussion). These general assumptions trans- 
late into \ogP{0\p) = Vf{0,p) and logP(p) = Vg{p,a) 
where / and g are large deviation functional whose de- 
tailed expression is not known in general, but are expected 
(generically) to be regular and have locally quadratic max- 
ima. We claim that these general results also hold for 
two-point correlation functions C(t, i-l-r), provided r is so 
small that p can be considered constant within the interval 
[t, t + r], such that the dynamics is slaved to a single value 
of p during this time interval. In this case, C{t, t + r) can 
be treated as a single time observable and we also have 
\ogP{C\p) = Vf{C,p). 

The dynamical susceptibility introduced in the con- 
text of supercooled liquids corresponds to the fluctuations 
5C = C - {C): Xi = N{SC^) where TV = {p)V and (•) 
stands for ensemble average. It can be easily computed 
from the large deviation functional by expanding around 
the most probable values, leading to Gaussian integrals. 
Following the discussion in [13], one obtains the following 
general relation: 



X4 = N {dC^)\+N{6p^) 



{ d{C{t,0)) 
V d{p) 



(1) 



This exact equation expresses the dynamical susceptibility 
as a sum of fluctuations in a restricted ensemble where the 
density is strictly flxed plus a term corresponding to the 
dynamical fluctuations induced by density fluctuations. 
Dropping the flrst term, which is clearly positive, one finds 
the inequality: 



X4> 



NjSp^) f d{C{t,0)) 

aiog(p) 



(2) 



In the following we shall measure directly all the quan- 
tities that appear in Eq. (2) and evaluate both the RHS 

and the LHS of the inequality in a horizontally vibrated 
amorphous assembly of bidisperse hard disks in order to 



examine whether the RHS provides not only a bound but 
a quantitative estimation of the LHS. 

Experimental results. — For tlu^ purpose of clarity, 
we flrst recall the essential features of our experimental 
system, which have been already described in [15]. A 1:1 
bidisperse monolayer of 8500 brass cylinders of diameters 
dsmail = 4 ± 0.01 mm and dbig = 5 ± 0.01 mm laid out on 
a horizontal glass plate. The plate is subjected to in-plane 
oscillation at a frequency of 10 Hz, with a peak-to-peak 
amplitude of 10 mm. The grains are conflned in a cell 
that is flxed in the laboratory frame. The cell has width 
L w 100 dsrnaii, ^ud its length can be adjusted by a lat- 
eral mobile wall controlled by a pm accuracy translation 
platen, which allows us to vary the packing fraction of the 
grains assembly by tiny amounts {S(j)/4) ^ 510"**). The 
stroboscopic motion of a set of 1500 grains in the center of 
the sample is tracked by a CCD camera synchronized with 
the plate. The protocol, detailed in [15], leads to reversible 
stationary states for all studied packing fractions, which 
range from (j) = 0.8402 to ^ = 0.8454. For each pack- 
ing fraction, the grains trajectories fi{t) are recorded for 
twenty minutes. The averages are performed over the cor- 
responding 10000 conflgurations of the packing. Lengths 
are measured in dg„iaU units and time in cycle units. 

At the jamming density, (pj ~ 0.842, where mechani- 
cal rigidity sets in, long range correlations develop both 
in time and space and give rise to a strong peak in 
the dynamical susceptibility [15]. Noticeably, the indi- 
vidual displacements themselves are much smaller than 
the grain size. A sketch of the behavior of the diffusion 
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is displayed in 



Fig. 1 (for a detailed discussion see [15]). Remarkably, 
the early subdiffusive regime connects to the late diffu- 
sive one through a transient super-diffusive regini{\ This 
regime, observed at characteristic time Tsd, emerges when 
approaching the jamming transition. A consequence of 
this striking behavior is that the variation of the diffusion 
length with the packing fraction, Sa^r) is larger at TgD 
leading to a non monotonic variation of (5(T(r) as a func- 
tion of time. The similarity with the behavior of xa is not 
accidental. We shall show in the following that the two 
turn out to be strictly related. 

We now concentrate on the evaluation of both terms of 
the inequality (2). For that purpose, we shall focus on a 
particular correlation function deflned as: 



(3) 



where a is the length scale over which the motion 
is probed. The tagged particle correlator Q (a, r) = 
{C{a,T)), when represented as a function of the probe 
a for a given r, is characterized by a decay length C,^ (r). 
This length scale is found to be proportional to the diffu- 
sion length (70 (t) , with the same constant of proportion- 
ality for all the investigated densities (Fig. (l)-Middle). 
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Fig. 2: Top: Shape functions X4 (D) and 
as a function of scaled variable x for = 
0.8402 (red) ,0.8417 (black) ,0.8426 (blue). The index n 
signals that the function has been normalized by its max- 
imum for comparison. Middle: Amplitudes (■) and 
dC,4, (•) as a function of lag r for three packing fractions 
= 0.8402, 0.8417, 0.8435. Bottom left: density fluctuation 
(Jp'^) (K) as a function of the radius R of the disk within which 
it is evaluated, for = 0.8417. The black line corresponds 
to the expected behavior. Bottom Right: Normalized 
density fluctuations N {5 p^) / {p)"^ as a function of packing 
fraction. The blue line is the interpolation we used to evaluate 
the bound. 



Fig. 1: Top: Sketch for the non-monotonic behavior of d(^^ . 
The diffusion length cr^, proportional to C,^, overshoots at 
intermediate times, giving rise to a temporal maximum in 
its density derivative. Middle: Q (a, r) as a function of 
scaled variable C0 (r) /a for r = 12 (O) , 108 (V) , 1047 (O) and 
= 0.8402 (red), 0.8417 (black), 0.8426 (blue). Lower insert: 
rough functions; upper insert: characteristic length scale C,^ 
vs. (J0 for all lags and densities. The dependence is linear and 
the slope C,^ — 1.62cr^ does not depend on 0. Bottom: Shape 
function X4 as a function of scaled variable (r) /a for the 
same set of data. Lower insert: X4 (^i ''') ^ ^ function of a for 
the three lags above; upper insert: the amplitude (r) as a 
function of lag r for the three densities above. 



The correlator Q hence follows the approximate scaling: 
Q{a,T)^Q{x) with a; = ° , (4) 



which accounts for a corresponding scaling of the displace- 
ments distribution probabilities. Perhaps surprisingly, the 
dynamical susceptibility X4 — N{6C'^) partly inherits this 
scaling property and reads 

X4{a,T) ^ h4,{T)xA{x), (5) 

where we normalize the maximum value of X4 , reached for 
X = 1, to be equal to 1 (see fig. 1-Bottom). The shape 
function X4 fully encodes the dependence on the probe 
lengthscale, whereas the amplitude h^{T) encodes the lag 
dependence of the heterogeneity of the dynamics for each 
packing fraction and reaches a maximum near Tgjj. 

Plugging the scaling properties (4,5) into inequality (2) 
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and using the chain rule yields 

(6) 

In this expression, one obtains a clear factorization of the 
probe dependence, contained in the shape functions, re- 
spectively: 

and the lag dependence in the amplitudes 

These two sets of functions can then be compared in- 
dependently. The shape terms represented on top of 
Fig. (2) were obtained by averaging over time. They ap- 
pear to match very neatly, and the overall factor is close to 
Xi ~ 5dQ . This factor is the squared maximum slope of 
Q and would be equal to 4 if the statistics of the displace- 
ments, in addition to obey the above scaling, were also 
gaussian. We can now turn to the comparison of the lag- 
dependent amplitudes and dC'^ represented on the mid- 
dle panel of Fig. (2) for three values of the packing fraction. 
These amplitudes have been normalized by their maxi- 
mum for proper comparison, and extend over the range 
of lags where derivation can be accurately achieved. Sur- 
prisingly, we observe that the response term is non- 
monotonic and that these amplitudes reach their maxi- 
mum for very similar lags. As stated above, the existence 
of the maxima in the response term is directly related 
to the existence of the superdiffusive regime: as readily 
visible from the sketch of Fig. 1 (Top), it confers a char- 
acteristic bell shape to da ^{t)/ dip cx dC^^{t) / dcj). 

The last step in investigating the quality of (2), is to 
compare the absolute ^-dependent bound 

= ^^^;^max,aC| max^Sg', (9) 

and the maximum dynamical susceptibility (recall that 
maxa;X4(a:) = 1): 

Xl {(f)) = max^/i0 (r) . (10) 

We have measured the density fluctuations (Sp^) within 
circular regions of increasing radius R and asymptotically 
found the expected behavior {Sp"^) {R) cx R~^ cx 
as can be seen on the insert of the bottom left panel of 
Fig. (2). As expected from a naive compressibility argu- 
ment, the corresponding intensive normalized amplitude 
N{5p^)/{p)'^ is a decreasing function of the packing frac- 
tion. The final evaluation of the terms (9) and (10) to 
be compared is shown on Fig. (3) for packing fractions al- 
lowing proper evaluation. The lower bound B is observed 
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Fig. 3: Final comparison of the maximum dynamical suscep- 
tibility Xa C) with the bound B (•) for the studied packing 
fractions. 

to match the dynamical susceptibility quite accurately, in 
particular it reproduces remarkably well the peak of dy- 
namical susceptibility across the jamming transition. The 
fact that it is found to be larger than Xi for high densi- 
ties is probably due to a loss of accuracy in our numerical 
derivatives and in the evaluation of the density fluctua- 
tions in that regime. 

Discussion. — Let us now discuss these results phys- 
ically. First, we have uncovered a strong scaling property 
in the wavelength dependence of both the dynamical sus- 
ceptibility and response, which emphasizes the fact that 
dynamical fluctuations at a certain time scale are maxi- 
mum when probed at the corresponding diffusive length 
scale. However, the remaining prefactor (t) still ex- 
hibits a maximum as a function of r. Therefore, evalu- 
ating this quantity at fixed probe-length or wave-vector 
does not in itself lead to the identiflcation of the dynamic 
fluctuation time. Second, we have found that the tem- 
poral behavior of the bound, driven by the derivative of 
the diffusion length with respect to the packing fraction, 
is a non-monotonic function of time. This behavior is due 
to superdiffusion, as sketched in Fig. 1. The fact that 
both sides of the inequality reach their temporal maxi- 
mum at comparable lags establishes a direct link between 
dynamical heterogeneity and anomalous autocorrelation 
of the grains displacements that leads to superdiffusion. 
The moving regions then appear to be disordered spa- 
tially but with persistent currents with typical life time 
TsD- Finally, the lower bound was found to account for 
a substantial fraction of the equality (1). This gives sup- 
port to the hypothesis we made to describe the statistics 
of the system, i.e. the absence of long-range correlations 
in the density field, and suggests the existence of a large 
deviation function for the intensive observables. Perhaps 
more importantly, the quality of the bound shows that 
density fluctuations is the main determinant of the statis- 
tical properties of the dynamics in our system, which indi- 
cates that the intensity of energy injection plays a minor 
role in this range of packing fractions. Finally, the bound 
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is found to be a very accurate indicator for the maximum 
amplitude reached by the dynamical susceptibility X4 at 
each density. 

Altogether the present work provides with the first 
experimental investigation of a dynamical fluctuation- 
response relation which gives a lower bound on the am- 
plitude of dynamical heterogeneities. The bound repro- 
duces the dependence of the dynamical susceptibility on 
probe-length, lag time, and correctly captures its density 
variation and amplitude. These important results clearly 
establish the relevance of this bound for characterizing 
dynamical heterogeneity in an athermal system. It would 
now be interesting to clarify the respective contributions of 
temperature and density in the context of thermal glasses. 
The next step would be to extract the "fixed density" con- 
tribution to Xi in Eq. (1) and understand in detail its 
physical origin. 
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